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o : 

CN| ■ To the first two orders, discrete holomorphic functions are whether complex 

O ' linear or Mobius transformations. We discretize this notion to complex functions of 

D . vertices <0o of an oriented quad-mesh <). First, fix a direct embedding of this quad- 

Ph ' mesh into the complex plane by a complex function of the vertices Z : <£>o —> C, 

y—i • corresponding to the identity map. We will say that another function is linear 

t-H ' holomorphic (with respect to Z) if and only if the ratio along the diagonals are 

the same as Z on each quadrilateral: 

Si , , f(y') - m Z(y') - Z(y) . 

^ V(x, y, x , y ) e 0,, /( ^ _ = ^ _ ^ = » P M ■ 

' Likewise we will say that a function is cross-ratio preserving if the cross-ratio on 

d . each quadrilateral is the same as the one given by Z: 

Si fix) - m fW - f(v') _ Z{x) - Zjy) Z(x') - Z(y') 

'— '; /(y) - /(x') /(t/') - /(*) Z(y) - Z(y') - Z(x) 'CW.w')- 

An important class of cross-ratio preserving maps are given by circle patterns 
, with prescribed intersection angles. In the form of a Hirota system, a connection 

between the two notions can be described: A discrete function F is cross-ratio 
, preserving if its exterior differential can be written, on each edge (x,y) G <C>i 

; as F(y) - F(x) = f (x) f (y) (Z (y) - Z{xj) =: J {Xty) f dZ for a function /. The 

£N| ■ constraint on / is the Morera equation: § f dZ = around every quadrilateral. 

| Looking at logarithmic derivatives f £ = fx (I + eg) of / that still form a Hi- 

rota system, one finds that g should be linear holomorphic with respect to F: 
. . , 9< ! y ,\~ 9 ^ v \ = tjtit— tft\- Linear constraints can as well be reformulated as a Mor- 

g(x')-g(x) F(x')-F(x) 

, £^ | era equation § f dZ = not for a multiplicative but for an additive coupling 

' between functions and I-forms: J^ x » / dZ := £M±£faj (Z(y) — Z{xj). 

^ ■ We will be interested in so called critical reference maps Z composed of rhombi, 

the length 5 = \Z(y) — Z(x)\ is constant for all edges (x, y) £ Oi- Then these two 
constraints are integrable, meaning that they give rise to a well-defined Backlund 
(or Darboux) transform: given a linear holomorphic rcsp. cross-ratio preserving 
map / and a starting point O £ v o, one can define a 2-paramcters family of 
deformations of / that still fulfill the same (linear resp. cross-ratio) condition and 
starting value [i at the point O. This is done by viewing the solution / lying on 
a ground level and building "vertically", above each edge (x,y), a quadrilateral 
on which the same kind of equation will be imposed: { Xai [ x } = f fo| 

H 1 h,t,{y)-f(x) \+Z{y)-Z{x) 

for the linear constraint, resp. ^Vw/^ / {v }~JY iv } \ = t n for the 

' 1 f(x)-f(v) f\,»{y)-f\,»(x) (z(y)-Z(x)) 

cross-ratio constraint. 

In [5], based on the point of view of integrable systems, we define discrete 
holomorphicity in Z d , for d > 1 finite, equipped with rapidities (ai)i<i<d and 
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the tools of integrable theory yields interesting results like a Lax pair governing 
a moving frame A) : Z d — > GL2(C)[A] and isomonodromic solutions like the 
Green function found by Kenyon [1] . 

An important tool of the linear theory is the existence of an explicit basis of 
discrete holomorphic functions. In the rhombic case, for a discrete holomorphic 
function /, the 1-form / dZ is holomorphic (it is closed and its ratios on two dual 
diagonals arc equal to the reference ratios) and can be integrated, yielding back 
a holomorphic function. Differential equations can be setup, producing explicit 
formulae for exponentials and polynomials which are shown to form a basis of 
discrete holomorphic functions. 

Lots of results of the continuous theory can be extended to this discrete settting: 
There exists a Hodge star: * : C k — > C 2 ~ k , defined by f, * a :— P( x ,x') L x x >\ a 'i 
the discrete Laplacian is written as usual A := dd* +d*d with d* : = — * d * , which 
reads as weighted differences around neighbours A f(x) = ^2 P(x,x k ) (f( x ) ~ f( x k))- 
The weights are given by the usual cotan formula. 

A wedge product defines an L 2 norm for functions and forms by (cx,/3) '■= 
//<>, a A * P- The norm of df is called the Dirichlet energy of the function /, 
E D \f) := \\df\\ 2 = {df, df) = | E(*,xOeA I P^ x ') \f( x ') ~ f( x )\ 2 ■ Thc informal 
energy of a map measures its conformality defect Ec(f) ■= \ \\df — i * df\\ 2 . They 
are related through E c {f) = E D (f) — 2A(f) just as in the continuous. 

The Hodge star decomposes forms into exact, coexact and harmonic ones, the 
harmonic being the orthogonal sum of holomorphic and anti-holormorphic ones. 
A Weyl's lemma and a Green's identity are found. 

Non closed 1-forms with prescribed diagonal ratios define meromorphic forms 
and the holonomy around a quadrilateral is called its residue. The compact case 
is covered with flat atlases of critical maps for a given euclidean metric with conic 
singularities. Meromorphic forms of prescribed holonomies and poles arc defined 
and are used to form a basis of thc space of holomorphic forms. It is 2<7-dimensional 
on a genus g surface, that is twice as large as the continuous case, defining two 
period matrices instead of one. This difference is explained by the doubling of 
degrees of freedom, and partially solved through continuous limit theorems: the 
two period matrices converge to the same limit when refinements of quad-meshes 
for a given Euclidean metric with conic singularities are taken. Every holomorphic 
function can be approximated by a converging sequence of discrete holomorphic 
functions on refinements of critical quad-meshes. 

The Green function and potential allow one to setup a Cauchy integral formula 
giving the value at a point (in fact its average at two neighbours x, y) as a contour 
integral: § dD f dG x , y = 2zir 

We define a derivation with respect to Z by d : C°(0) - ¥ C 2 ((}) with df = 
Ut a t' n'\ i-» i (L friz] — (■f( x ')-f( x ))(.y'-y)-( s '- s )(f(y')-f(.y)) 

\)f L iti ) " 2A(x,y,x',y') J J UZj I ~ (x' -x)(y' -y)-(x' -x)(y' -y) > 

(x,y,x\y') 

where Z(x) is simply written x; and likewise df with fdZ. A holomorphic func- 
tion / verifies df = and df{x,y,x',y') = — = £^2J^£M. The Jacobian 

y y xx 
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J = |<9/| 2 - \df\ 2 relates the areas // df A df = J JJ dZhdZ. 

(x,y,x',y') (x,y,x',y') 

Following Colin de Verdiere and Kenyon, a geometrical 
interpretation of linear discrete holomorphicity is enlight- 
ening. As circle patterns with prescribed angles can be 
checked by eye, so can be a linear holomorphic map: The 
quad-mesh <>, when bipartite, decomposes into two dual 
graphs T and T* whose edges are dual diagonals of each 
quadrilateral. Around each vertex x £ To, there is a poly- 
gon, image of the dual face x* s T% by the reference map 
Z. Consider the identity map as a picture of all these 
polygons shrunk by a factor half. It represents both dual ~^^Sr^XYX^ 
graphs at the same time as matching polygons. A map ^ ►iT^ 

f : <^>o — > C is discrete holomorphic if and only if every 

polygon x* , centered at f{x), scaled and turned accord- -phe discrete exponential 
ing to df(x), form into a polygonal pattern of the same as a polygonal pattern on 
combinatorics as the reference polygonal pattern, made the triangular/hexagonal 
of similar polygons. lattige 

The dilatation coefficient of a discrete map / is defined as Df := ^ • 
will call / quasi- conformal when Df > 1, that is < |<9/|. It can be written in 
term of the comnler dilatation ■ n t — ^ - 

term or tne complex auatauon. fi f - 8f - {f { x ')-f{ x )){y'-y) - (x'-s) (f(y')-f(y)) ■ 
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